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Method and unit for substracting quantization noise from a PCM signal ^ ^ 

2 5. 01. 2CG2 
@ 

The invention relates to a method and a unit for substracting quantization 

noise from a PCM signal. 

Pulse code modulation (PCM) is a widely used technique to represent e.g. 
audio signals in the digital domain. Linear PCM uses evenly distributed quantization levels 
whereas the logarithmic PCM uses a distribution where small-signal amplitudes are 
quantized using more closely spaced levels than large-signal amplitudes. An example to 
linear PCM is the CD-audio, which uses 16 bits to represent uniformly distributed signal 
levels. 

Fig. 5 illustrates the quantization of a linear coded PCM signal amplitudes 
(curve (a)) and of two logarithmicaUy coded PCM signal amplitudes (curve (b) and (c)). 
Cuve b illustrates the quantization of coded PCM signal ampUtudes, wherein the jx-law 
has an input-output magnitude characteristic of tixe form: 



where l^j is the magnitude of the input. |yl is the magnitude of the output and// is aparameter 
that is selected to give the desired compression characteristic. 

In Fig. 5 curve c) illustrates the quantization of A-laW coded PCM signal 
ampUtudes. wherein the A-law has an input-output magnitude characteristic of the form: 



lH4og^ (2) 
l+logA 

wherein the meaning of y and s is as mentioned above and A is constant The prior art jx-law 
and A-law coders are standardized by the Comit6 Consultatif International de T616f onique et 
T616graphique CCHT as the G.711 64kbit/s PCM. The quantizer characteristics of tiie G.71 1 
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are disclosed in: N.S. Jayant and P. Noll, "Digital Coding of Waveforms", Englewood Cliffs 
NJ: Prentice-Hall. 1984. 

The quantized sample values of a signal may be modeled mathematically as 

5 s^[n\ = s[n\-\-q[n], (3) 

wherein S*[n] represents the quantized value of an original unquantized sample value s[n] 
and q[n] represents the quantization error resulting from said quantization. The quantization 
error can be interpreted as undesired additive noise decreasing the quality of the original 
10 sample values s[n]. 

Thus, it is the object of the invention to provide a method and a unit for 
subsfcracting said quantization noise from a PCM signal. 

This object is solved by a method comprising the steps of : 
calculating for each frame of said PCM signal a quantization noise level Bq according to the 
15 following equation: 



20 



^ ^' ;s 12 



wherein 

n : indicates a specific sample of the PCM signal; 

S*min[n] : represents the minimum quantization noise level for a specific sample value 

s*[n] of said PCM signal; 
25 S*max[n] : represents the maximum quantization noise level for the specific sample 

value s*[n] of the PCM signal; 
w[n] : represents a window-function; and 

W : represents the number of samples per window; 



30 and 



substracting the quantization noise as represented by said quantization noise level Bq from 
said PCM signal. 
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It is the essential idea of the invention to calculate said quantization noise 
levels Bq for each frame according to the claimed formula and to use any suitable 
(background) noise substracting method known in the art - perhaps after carrying out sUght 
modifications - to substract the quantization noise from the PCM signal. 
5 Throughout the whole description * indicates the content of quantization noise 

in the co-ordinated quantity or signal. 

The window-function may be chosen arbitrarily. 

The claimed formula requires the minimum quantization level S*fflin[n] as 
well as on the maximum quantization level S*^[n] of the sample values s*[ii] of the PCM 
10 signal. 

According to a first advantageous embodiment of the invention the minimum 
quantization level S*«ia[n] as well as the maximum quantization level S*„««[n] are known 
from the quantization method applied for generating the PCM signal. 

According to a second embodiment said quantization levels S^^nM and 
15 S*«a«[n] are not known and must therefore be predicted according to the formulas claimed in 
claim 3. 

Advantageously neither in the first nor in the second embodiment any 
additional information must transmitted together with the PCM sample values s*[n\ to enable 
the conduction of the claimed method. Hiis is because in the first embodiment the 
20 quantization levels and in the second embodiment reference levels i can be assumed to be 

known (see the last two paragraphs above) when the method for substracting the quantization 
noise acconiing to the invention is carried out Because of that advantageously bandwidth is 
saved. 

Preferably, the noise substraction is based on a background noise substracting 
25 method and system as described in the non-prepubUshed European Patent appUcation having 
the application number BP 01201304.1 of the same appUcant. The appUcation of said system 
is described later as preferred embodiment. The whole content of said patent appUcation shaU 
be regarded as part of the disclosure of the present appUcation. 

Fmally. the above-identified object is solved by a noise substi-acting unit 
30 according to the subject matter of claim 9. The advantages of said unit correspond to the 
advantages mentioned above referring to the claimed method. 

Preferably, the claimed unit and metiiod are appUed at tiie decoder side, i.e. 
just before tiie PCM signal is decoded again. However, nevertheless it might also be appUed 
at tiie encoder side. 
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The description is accompanied by five figures, wherein 
Rg. 1 shows a quantization noise substracting unit according to the invention 

including the background noise subtracting unit; 
Fig. 2 shows an example for illustrating the correlation between reference values I 

and sample values S*[n]; 
Fig. 3 shoes a detailed illustration of a background noise substracting unit; 

Fig. 4. 1 . shows the spectrum of a clean male voice; 
Fig. 4.2. shows the |i-law coded version of the spectixim shown in Fig. 4.1.; 
Fig. 4.3. shows the ^i-law coded version of the spectirum shown in Fig. 4.1. after 

quantization noise substraction according to the invention; 
Fig. 4.4. shows the usage array u[i] generated for applying the quantization noise 

substraction mediod to the spectirum shown in Fig. 4.2; 
Fig. 5 showing quantization levels for Hnearly and non-linearly coded PCM-signals 

as known in the art. 

In the following a preferred embodiment of the invention will be described by 

referring to Figures 1-4. 

Fig. 1 shows a quantization noise substracting unit 100-700 for substracting 
quantization noise from an input linearly or logarithmically coded PCM signal. 

First of all said unit comprises a quantization noise level calculating unit 100 
for calculating a quantization noise level Bq for each firame of said PCM signal under 
consideration of both a minimum and a maximum quantization noise level S*mm[n], S*n„«[n]. 
More specifically, the quantization noise level Bq is constant for each frame of said PCM 
signal and calculated individually for a specific frame according to die following equation: 

A| /i=o 12 



B 



wherein 

n : indicates a specific sample of the PCM signal; 

S*min[n] : represents the minimum quantization noise level for a specific sample value 

s*[n] of said PCM signal; 
S*max[n] : represents the maximum quantization noise level for tiie specific sample 
value s*[n] of die PCM signal; 
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v,rn] ' represents a window-function; and 

W : represents the number of samples per window comprising at least the specific 

frame for which Bq is currently calculated. 

{(.9* .rnl-S*-[nD-w[n]y_^ represents tben^an 

In equation 4 the term ^2 

square value of the quonfeaaon em,M.e. fl» energy of the quantizaaon noiso in to current 
frame. 

The frame size is typically 2ms. 

According to a first embodiment said quantization noise levels S*^„[n], 

S*maxtnl are assumed to be known. 

However, according to a second embodiment of the invention these 
quantization levels S*.„[n], S*^[n] are not assumed to be known, e.g. because the ong^nal 
signals might havebeen quantized with different quantizers and there may have been ex^^^ 

scalings appHed on the samples before or after coding. In these cases the quantization level. 

have to be predicted. , . ^ ' 

To predict the quantization levels, the entire original signal, heremafter also ■ 
xeferredtoaswaveform.isscam.edtoreceivethesamplevaluess*[n] of the quantized 
waveform. Subsequentiy.acheck is made, if all possible representation levelsipredefmed 

due the quantizationmethodappUedto the waveform arerepresented(used)byatle^^^^^^^^^ 
saidsamplevaluess*[nl.Foral6-bitlinearPCMquantizationmethodthereare2 =65536 
possible representation levelsiand the result of saidcheckis described inaso called usagp 

table u[i] with 65536 values. T^e values of said usage table u[i] are either 1 (if there is at 
least onesamplevalueof the whole waveform corresponding to the current representation 

level i) or 0 (if there is no sample value of the whole waveform corresponding to the current 
representation level i). Expressed mathematically 



15 



20 



25 



30 



= Otherwise) 



.2'''<i<2'-' (5) 
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In formula (5) L is the wavefoim length, i.e. the number of samples of the 
whole waveform, and N is the number of bits used by the Unear PCM quantization; in the 
above example N=16. This operation can be simplified by initializing all values of u[i] with 0 
and using formula (6): 

M[5*[n]] = l, 0<n<L (6) 



10 



15 



Now equations (5) or (6) are used to find special representation levels 
and w^n; required for predicting the desired quantization levels. More specifically. i„dn[n] 
is defined as that representation level which is - startet from i[n] - the next smaller non-zero 
representation level for which u[n]=l and ir^[n] is defined as that representation level which 
is - started from i[n] - the next bigger non-zero representation level for which u[n]=l. 

With the help of said co-ordinated representation levels the desired 
quantization levels are predicted or calculated according to the following equations: 



20 



25 



wherein 
i 

i[n] 



S*niin = iM - (iM - i«in[n] ) / 2 
S'nm = i[n] + (wW - i[n] ) / 2 



(7) 
(8) 



: represents one out of a pluraUty of predefined possible representation levels; 
: represents that predefined representation level which corresponds to tiie 

sample value s*[n]. 
Fig. 2 illustrates the calculation of the usage array, of said special 
representation levels i„anfnj, inuJn] and of the desired quantization levels S*«i„[n], S*™ax[n] 
for a given waveform consisting of 1^9 sample values S*[n] which are in tiie example 
distributed as shown in Tab. 1: 



30 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


s*[n] 
=i[n] 


3 


2 


-1 


-1 


2 


1 


1 


3 


-3 



Tab. 1 
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In Fig. 2 there are provided 8 different possible representation levels i= -4, 

10, 1 3 resulting from a N=3 bit Unear PCM quantization. The usage array u[il calculated 

according to equation (5) or (6) for the example of Tab. 1 is: 

5 

u[-3] = u[-l] = u[l] = u[2] = u[3] = 1; and 
u[-43=u[-2] = u[0] = 0 

Now, from said usage array the reference levels iminW and inaxM are 
10 calculated e.g. for n=2 with s*[2]=i[2] =1 - according to the definition given above - to: 

imax[2]= 2 and ij^[2] - -1 

and the desired quantization levels result -according to equations (7) and (8) to: 

15 

s*„^=i-(i-(-i)y2=o 

S'fflax=l+(2 -1)/2 = 1,5 

Special care has to be taken when i[n] co-ordinated to s*[n] is the first or last 
20 non-zero value in the usage-array u. One option is to allocate two extra array values ul-Z"''- 
1 J and ufl"-' J, set them both to 1 and always continue the search until a non-zero value is 
found. 

Another question is what to do when some valid quantization levels aren't 
used in the waveform. Especially waveforms shorter than 1 second will result in a usage- 

25 anay with many unused (but valid) quantization levels. Since this usually wiU give a wrong 
result for the quantization step-size at boundaries (when the waveform is at its maximum or 
minimum value), the effect on the measured level of quantization noise will be small. 
Although an algorithm could be used to predict the unused vaUd quantization levels, the 
added value of the extra complexity is not worth the performance gain. In practice, it seems 

30 to be enough to assume that the used signal levels in the waveform are the only valid levels. 

After the quantization noise level Bq has been calculated either with the 
known or with the predicted quantization levels s*^n and s*^ it is input into a background 
noise substracting unit BNS 200 which is also part of said quantization noise substracting 
unit. 
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Fig. 3 shows a detailed iUustration of said background noise substracting unit 
200 which is embodied to substract the quantization noise from a PCM signal in the 
frequency domain. It comprises a magnitude forming unit 210 for forming the magnitude 
1S*M1 of the spectrum S*[k] of the PCM signal, a signal-to-noise ratio calculating unit 220 
for calculating the signal-to-noise ratio of the PCM signal according to the formula: 

SNR[k] = lS*[k]imq. C9) 

wherein Bq is noise level representing the level of quantization noise in said PCM signal. 

Said signal-to-noise ratio SNRM is input into a filter update unit 230 for 
calculating a filter magnitude F[k] according to a predefined filter algorithm based on at least 
one filter update parameter. Said filter magnitude F[k] is input into a calculation unit 240 for 
calculating an output spectrum S^M by multiplying both the real part R{S*[k]} and tiie 
imaginary part I{S*[k]} of the spectmm S*[k] with said filter magnitude F{k}. The filter 
update parameter and thus also the filter magnitude F[k] are preferably adjusted such that 
output spectrum S^'M is at least substantially free of quantization noise. 

Because said background noise substracting unit 200 is embodied to carry out 
the substraction of the quantization noise in the ftequency domain, the quantization noise 
substracting unit further comprises several units for computing the spectrum S*[k] of the 
PCM signal to be input into said background noise substracting unit 200 and for transforming 
the spectrum S^M output by said unit 200 back into a time domain signal. These units are 
described in the following paragraphs. 

Referring back to Fig. 1 the quantization noise substracting unit further 
comprises a first window block 300 for generating a weighted PCM signal s*w[n] by 
weighting the firames of the input PCM sgnal with a first window wl[n]. Said weighted PCM 
signal s*w[n] is input into a spectmm computing unit 400 for computing the spectrum S*[k] 
from said weighted signal. Preferably, the spectnam is computed by using the Fast Fourier 
Transfomiation algorithm. The spectrum S*[k] is than input into the background noise 
substracting unit 200 and processed therein as described above by referring to Fig. 3. 

' The spectrum S'^M output by said background substracting unit 200 is input 
into a le-transformmg unit 500 for being transformed back into a signal s^[n] in the time 
domain. In tfie case that the spectrum S*[k] has been calculated by using FFT the re- 
transformation is done by using the inverse FFT. 
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The quantization noise substracting unit may (uither comprise a second 
windo* block 600 for generating a wd^tcd output signal s'.M by weighting tlte signal 
sVl received after the transformation with a second window W2M and an overlap and add 
»it700forindividuallycalculadngafinal output signal Slltl] for each frame of U« 
, PCM-signal ftom said wei^ted output signal s'„[n] such that the trandtion between two 
successive output frames is smoothed. The first and the second window Wl and ^Vz are 
preferably idenUcal. Ttoy are independent of the window-function wM used for calculattng 
d,e qoantizadon noise level Bq. Hov,ever, it is preferred to ose the same function tor w W ^ 

10 Thepiovision of the overlap and add unit 700 only improves the quality of the 

outputsignalbutitisnotabsolutelynecessaryforcarryingoutthequantizationnoise 
substraction; consequently . it xnight be left out. The san. applies to *e first and 
window block 300 and 600 v.hich are both only provided to improve the operation said 
overlap and add unit 700. - 

15 Fig. 4 shows an example illustrating the good signal qoaUty achieved when 

applying the claimed method to a Maw coded speech signal. i 

Mote specifically. Pig. 4.1 shows the spectrum of a clean, i.e. ftee of 
quantizationnoiscmalevoicerecordingwith the utt.r«>ce: "Any dictionary will ^ve at least 
twenty different deflnitions for it". Hg. 4.2. shows the ^^law coded version of the spectrum 
20 shown in Hg. 4.1. Prom that Hg. 4.2. it can be seen that said ^aw coding causes 

quandzatton noise added to the clean ap^tmm of Hg. 4.1. In Kg. 4.2. the qnantizaUon no.se 
is represented by the grey-shaded areas 410; the corresponding areas m Hg. 4.1 are clean. 
Fig 4 3 shows the spectrum from Hg. 4.2 after the quantization noise substracUng method 
,^g,„ the invention has be^apphed. I. can been seen that MS spectrum substantially 

25 corresponds to original clean spectrum shown in Hg. 4. 1; that means that the method 
according to the invention is well able to clean the spectrum of a PCM signal from 
quantization noise. Hg, 4.4 shews the usage array n[i] as generated by the method accordmg 
to the ptesent invention for reducing quantization noise from the spectrum shown in Hg. 4.2. 
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^PO - DG 1 

CLAIMS: 

25. 01 2GQ2 



1. A method for substracting quantization noise from a pulse code modulated 

PCM signal being segmented into frames, comprising the steps of: 
calculating for each frame of said PCM signal a constant quantization noise level Bq 
according to the following equation: 



" 12 



10 wherein 



n : indicates a specific sample of the PCM signal; 

S*nBn[n] : represents the minimum quantization noise level for a specific sample value 

s*[n] of said PCM signal; 
15 S*Bax[n] : presents the maximum quantization noise level for the specific sample value 

s*[n] of the PCM signal; 
w[n] : represents a window-function; and 

W ' represents the number of samples par window; 

20 and 

substracting the quantization noise as represented by said quantization noise level Bq from 
said PCM signal. 

2. The method according to claim 1 , characterized in that the minimum 
25 quantization level S*mm as well as the maximum quantization level S*max are known. 

3. The method according to claim 1, characterized in that the mimmum 
quantization level S*^*, and the maximum quantization level S*„ax are predicted according to 
the following equations: 
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S*roin=i[n]-(i[n]-ix„in[n])/2 
S=max = + (imaxN - i[n] ) / 2 



wherein 
i 

i[nl 

iminLn] 



: represents one out of a plurality of possible representation levels predefined 
due to the specific PCM quantization method applied to an original signal; 
: represents that predefined representation level which corresponds to the 
sample value s*[n] for a specific n; 

: represents that representation level which is - startet from i[n] - the next 
smaller non- zero representation level for which u[nl=l; 
: represents that representation level which is - startet from i[n] - the next 
bigger non- zero representation level for which u[n]=l; 



with the usage array u[i] being defined to: 



r 



u\i\ = min 



n=0 



0. 
1, 



5 * [ra] 5!: i 

otherwise 



wherein 

L : represents the number samples of the whole PCM-signal; and 

N : represents the number of bits used for quantizing an original sample value 

by using PCM to generate the PCM sample values s*[n]. 

4. The method according to claim 1, characterized in that the substracting of the 

quantization noise represented by said quantization noise level Bq from the PCM-signal is 
carried out in the frequency domain according to the following steps: 

computing the spectrum S*[k] of the PCM signal s*[n] and forming the 

magnitude \S*\k]\ thereof; 

computing a signal-to-noise ratio SNR[k] of said spectrum S*Ilc] according to: 

SNR[k] = iS'^Mj/Bq; 
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calculating from said signal-to-tioise ratio SNR[k] a filter magnitude F[kl 
according to a predefined filter algorithm based on at least one filter update parameter; 

calculating an output spectrum S^M at least substantially free of said 
quantization noise by multiplying both the real part R{S*[k]} and the imaginary part 
5 I{S*[k] } of the spectrum S*[k] with said filter magnitude F{k}; and 

transforming the output spectrum S\k} back into a signal s^M in the time 

domain. 

5. The method according to claim 4, characterized in that the filter update 

10 parameter and thus the filter magnitude FM are adjusted such that the quantization noise in 
the remaining output spectrum s''[k] is as low as possible. 

6. The method according to claim 4, characterized in tiiat it further comprises the 
steps of : 

15 . weightingtheframesoftheinputPCMsignal withafirstwindowwlM and 

calculating the spectrum S*[k] from said weighted signal; 

generating a weighted output signal s\[n] by weighting the signal s^En] 
received after the re-transformation with a second window W2[n]; and 

calculating a final output signal [u] for a current fi^e of the PCM-signal 
20 from said weighted output signal s^[n] such that tiie tiransition between two successive 
ouq>ut firames and is smoothed. 

7. The method according to claim 4, characterized in tiiat 

the computation of the specti^im S* [k] of the PCM signal is done by using a Fast Fourier 

25 Transformation FFT; and 

the re-transforming the output spectrum S^M back into a time domain signal sV] is done by 

using an inverse FFT. 

8. The metiiod according to claim 6, characterized in that the first and the second 
30 window wl and w2 are identical. 

9. A quantization noise subshticting unit for substracting quantization noise from 
a pulse code modulated PCM signal being segmented into fiames, comprising: 
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a 

si: 



quantization noise level calculating unit (100) for calculating for each frame of said PCM 
gnal a constant quantization noise level Bq according to the foUowing equation: 



B = 



11=0 12 



wherein 
n 



: indicates a specific sample of the PCM signal; 
S*mia[n] : represents the minimum quantization noise level for a specific sample value 

s*[n] of said PCM signal; 
S*m«[n] : represents the maximum quantization noise level for the specific sample 
value s*[n] of the PCM signal; 
: represents a window-function; and 
W : represents the number of samples per window; 



as 



and 

a background noise substracting unit (200) for substracting the quantization noise 
represented by said quantization noise level Bq firom said PCM signal. 

10. The noise substiracting unit according to claim 9, characterized in that it is 

located at a decoder's side. 
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The invention relates to method and unit for substracting quantization noise 
from a pulse code modulated PCM signal being segmented into frames. For achieving this it 
is proposed to first calculate for each frame of said PCM signal a quantization noise level Bq 
according to the following equation: 



n=0 12 



wherein n indicates a specific sample of the PCM signal, S*„*.[nl represents the minimmn 
quantization noise level for a specific sample value s*[n] of said PCM signal. S*„K«[n] ; 
represents the maximum quantization noise level for the specific sample value s^n] of th^ 
PCM signal, w[n] represents a window-function and W represents the number of samples per 
window. Subsequenfly, the quantization noise as represented by said quantization noise leyel 
Bq has to be substracted from said PCM signal, preferably with the help of a suitable 
background noise substracting system. 
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